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Abstract. We compare entanglement with the violation of a Bell inequalitiy using 
a geometric structure in the state space of bipartite qudits. Determining whether a 
given quantum state obeys or violates a Bell inequality is a high-dimensional nonlinear 
constrained optimization problem for which in most cases there is no analytic solution 
available. Based on an advantageous parameterization we are able to study numerically 
a special set of states of bipartite qudit systems, namely convex combinations of 
generalized Bell states defining a simplex. This simplex is ideal for studying properties 
of density matrices due to the fact that any element can be descriptively represented 
by a real vector. We find that the shape of the boundaries of separability and states 
violating the CGLMP-Bell inequality are significantly different. Moreover, we confront 
the strength of violation with a recently proposed entanglement measure. We show 
that the implied orders for these density matrices are different, i.e. in general there 
exists no monotonous relation between a Bell inequality violation and the amount of 
entanglement. 
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1. Introduction 

The fact that quantum physics contradicts local-realism is one of the most seminal 
discoveries. In his pioneering work Ref. [I], John S. Bell showed that the statistical 
behavior of a bipartite qubit system is irreconcilable with any local-realistic theory when 
the system is in a singlet state. He revealed that the correlation probabilities of such 
theories satisfy certain constraints which can be violated if the setup is handled quantum 
physically. Since then much effort has been made to fully understand the origins and 
conditions that allow to demonstrate and experimentally test this contradiction known 
as quantum nonlocality. A variety of so-called Bell inequalities are nowadays known. 
Recent Bell inequalities also cover quantum systems with more than two degrees of 
freedom, e.g. Ref. [2] and/or more than two parties Refs. [31 HJ El E]- 

Yet, still little is known about the relation between entanglement and nonlocality 
apart from the fact that any violation can only be achieved by entangled states. It is 
not understood why there are entangled states that do not violate any currently known 
Bell inequality. The existence of such states has raised the question of how far these 
inequalities can be improved and if for a given system it is possible to specify the most 
powerful inequality or a set of those. In the end, are there entangled states allowing a 
local-realistic description even if LOCC (local operators and classical communication) 
and POVMs (positive operator valued measurements) are used? Moreover, a connection 
between the violation of a Bell inequality and the security of quantum cryptography has 
been shown, attesting to the importance of further investigations [7]. 

The conjecture that nonlocality and entanglement might not exactly be the same 
has been strengthened by the analysis of the CGLMP-Bell inequality [2J . In particular, 
it was shown that the maximal violation is achieved with non-maximally entangled 
states [8]. A similar feature was also found in Ref. [S] considering neutral K-mesons, 
which are two level systems oscillating and decaying in time. For this system entangled 
in strangeness though having only two degrees of freedom one finds no violation of the 
CHSH-Bell inequality if the system is prepared in a the spin singlet state copiously 
produced in accelerator experiments [10]. However, the maximal violation can be 
achieved if a non-maximally entangled state is used. It should be noted that the 
nonlocality of this particular system in high energy physics is considerably different 
from an ordinary bipartite qubit system and it has been shown for a certain choice of 
observables a relation between a Bell inequality and the tiny violation of a symmetry in 
particle physics, the CV violation (C. . . charge conjugation, V . . . parity), can be found, 
Refs. [HIE]. Another hint that entanglement and nonlocality are different resources is 
given in Ref. [13]. Herein, it is shown that the simulation of non-maximally entangled 
states via so-called non-local machines requires more resources than the simulation of a 
Bell state. 

In the present work we aim to study nonlocality and entanglement in the context of 
state space geometry. This concept is based on the idea of finding descriptive real vector 
representations for density matrices. In case of a single qubit such a representation is 
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given by the well-known three dimensional Bloch-vector with the state space given by 
the Bloch-sphere. Because of the high dimensionality of multipartite qudit systems 
one usually tries to find representations where local properties and correlations can be 
separated in different quantities. In all cases, when it comes to studying entanglement 
and nonlocality of a certain system it suffices to investigate representatives of equivalence 
classes of states which are equivalent with respect to local unitary transformations. In 
that way, for bipartite qubits is was found that correlations can be expressed with a three 
dimensional vector lying within a regular tetrahedron [HJ E]. For the investigation 
of bipartite qudit systems we utilize a proposed generalization of this tetrahedron, a 
simplex spanned by mutually orthogonal maximally entangled states [T6| \T7\ 118] . The 
goal is to visualize and compare the geometries of the boundaries of entanglement and 
nonlocality of low-dimensional subsections of this set of states. 

In order to do so, we choose certain classes of states for which via optimal 
entanglement witnesses we are able to determine the boundary between separable and 
entangled states. For the same classes we then determine all states which violate the 
CGLMP-Bell inequality. Determining, for a given state, whether it violates a certain 
Bell inequality or not is a high-dimensional nonlinear optimization problem which is a 
nontrivial task even for low dimensional quantum systems. The problem can at least be 
solved numerically when it is brought into an adequate form. In Sec. [3] we achieve this 
by utilizing the generalized Euler angle parameterization of the special unitary group 
SU(N) and removing redundant parameters. 

A visualization and comparison of the sets of nonlocal and entangled states for 
bipartite qubits is given in Sec. and for bipartite qutrits in Sec. [3 Finally in Sec. M 
we investigate the relation between the strength of CGLMP-Bell inequality violation 
and an entanglement measure, recently introduced in Ref. |19j . The results reveal that 
there exists no monotone function between both quantities. 

2. Foundations of Bell inequalities 

We begin by reviewing the foundations of general Bell inequalities. Any local-realistic 
description of a bipartite system with d outcomes and n measurement apparatuses on 
each side implies the existence of d 2n probabilities 



= j,...,A n = k;B 1 =l 



,...,B n = m)>0 



(1) 



with normalization 




P{A 1 = j, . . . , A, = k; B x = I, . . . , B n = m) = 1 



(2) 



j,.. .,m=0 



determining the statistics of the observables A±, . . . , A n , B±, ... , B n which can take on 
the values j, . . . , k,l, . . . ,m G {0, . . . , d — 1} [201 122]- As it is known, probability 
distributions of this form cannot reproduce the statistical behavior of composite 
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quantum systems 7i AB = C d <S> <C d when the system is in a certain state and the 
observables are chosen appropriately. In order to reveal this, one has to construct 
quantities / that consist of linear combinations of joint probabilities 

P(A a = x;B b = y) = 

d-i 

2j P(Ai = j, A a = x, A n = k; Bi = /,..., B b = y, B n = m) 

j,...,m=0 (except x,y) 

and show that they are bounded, i.e. 

I < Q (3) 

if ([1]) and (j2J) are assumed. If this bound can be exceeded when the system is treated 
quantum physically then the established inequality (jHJ) is called a Bell inequality. 

3. Determining nonlocal quantum states 

Generally, states that cannot be described by a local-realistic theory are said to be 
nonlocal. One way of detecting nonlocality of a certain state p is to show that it 
violates a certain Bell inequality. The set of nonlocal states is complementary to the set 
of local states whose elements allow probability distributions of the form ([TJ with (J5J) 
for any number of apparatuses n. It should be noted that locality of a state is not easily 
certified because, in principle, it has to be shown that any Bell inequality is satisfied or 
an explicit local realistic description has to be found. 

Our goal is now to choose a particular system and Bell inequality and to find a 
method to derive all entangled states which violate this particular Bell inequality. This 
is in principle a nontrivial problem. 

In detail, to show that a certain state p possesses nonlocal correlations one has to 
consider quantum mechanical joint probabilities 

P QM (A a = x;B b = y)= Tr(\x) Aa (x\ Aa ® \y) Bb (y\ Bb • p) (4) 

where {\x) A } and {\y) B } are orthonormal eigenvectors of the corresponding observables 
of Alice {A\, . . . , A n } and of Bob {B\, . . . , B n }- These probabilities are then inserted 
into the quantity /, which can then be rewritten in the form 

I = Tr(B lP ) (5) 

wherein Bi is called Bell operator. 

In order to examine if a Bell inequality is preserved or violated for a given p one 
has to be able to compute the global maximum of / with respect to all possible Bell 
operators 



max / = max Tr (Br 



(6) 
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This is in general a high-dimensional nonlinear constrained optimization problem for 
which analytic solutions are only known for a few special cases [23, 24J. In most cases 
as in the following investigations the quantity I has to be maximized numerically for 
each given p. 

Numerical methods can best be applied when an adequate parameterization with 
respect to the constraints is chosen and redundant variables are removed. In the 
following we show how this can be achieved: For all observables A a and there are 
corresponding orthonormal bases {|&)^ a } and {\k) Bb } determining the joint probabilities 

(A a = x,Bf, = y) and hence affecting the Bell operator. Any basis {\k)} 
of C d can be transformed into any other basis {!&')} of C d by means of a certain 
unitary transformation = {U \k)}. Consequently, one can choose an orthonormal 

computational basis {\k)} and seek the unitary transformations {UA a } and {UB b } such 
that I becomes maximal with {|^)^ a } = {UA a \k)} and = {Ub b \k)}- The 

transformations are elements of the unitary group U(d). Due to the fact that global 
phases do not affect probabilities P® M (A a = x,Bb = y) and that U(d) can be written 
as a semidirect product U(d) = SU(d) x £7(1), it suffices to regard the special unitary 
group SU(d). Without loss of generality one can use the generalized Euler angle 
parametrization [25l [26] of this group, wherein the Euler angles embody all degrees of 
freedom. For this kind of parameterization the diagonal generalized Gell-Mann matrices 



(k - l)k 



'k~2 



0'l-(*-i)l*-i>(*-i| 



(7) 



with 2 < k < d, and antisymmetric generalized Gell-Mann matrices 



A 



-z|0) (k\+i\k) (0| 



(8) 



with 1 < k < d— 1 are needed. In this notation the explicit form of the parameterization 
of an Ud G SU(d) by generalized Euler angles {c^} reads 



Ua 



'd-2 /d-x 



n (HA(k,j(x)) 



x=0 \k=2 



\ J exp(i\ n 2_ 1 a d 2_ {d+1 ^ n) ) 



n=2 



(9) 



with 



and 



M k , = exp(i\ 3 a(2k~3)+j(x)) ■ exp(i\ {k _ 1) 2 +1 a 2 (k-i)+j(m)), 



x = 

x-l 

2^2(d-x + l) x>0. 

1=0 



| Note that our notation differs slightly. We have chosen an indexing where all indices are increasing 



and where the sequence of the product is YiiLi A% — A\ ■ A2 ■ ■ ■ A 



N ■ 
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For the computational basisvectors \k) the last part of the parameterization only 
produces redundant global phases, i.e. Ylt=2 ex P(^n 2 -i a d 2 -{d+i-n)) \k) = e l ^ k \k), and 
can therefore be removed in our case. The numerical search can be restricted to the 
transformations 



d-2 /d-x \ 

U d = H(l[A(k,j(x))\ (10) 

x=0 \k=2 / 

For instance, for d — 2, 3, 4 one finds explicitly 
U 2 =exp(i\ 3 a 1 ) ■ exp(i\ 2 a 2 ) , 

U3 =exp(iA3a;i) ■ exp(i\ 2 a2) ■ exp(iA3Q;3) • exp(iA5a4) • exp^Asas) ■ exp(iA2«6) , 
U4 =exp(iA3«i) • exp(iA2«2) • exp(iA3«3) • exp(iAs«4) • exp(iA3a5) • exp(iAioa;6) 
•exp(iA3Q;7) ■ exp(i\ 2 as) exp(iA3ag) ■ exp(iAsaio) • exp(iA3an) ■ exp(?A2ai2) • 

The parameterization has a concise form, guarantees compliance with the constraints 
and reduces the number of involved parameters to d 2 — d for each observable without 
discarding any solution. Numerical optimization algorithms such as Differential 
Evolution [27], Simulated Annealing [28j or the Nelder-Mead method [2jJ] can now be 
used to maximize Tr(pBi) with respect to the remaining d 2 — d Euler angles «j. 

4. The CGLMP-Bell inequality 

We focus now on the CGLMP-Bell inequality [2J: If one considers the scenario of a 
bipartite system with each d degrees of freedom where the two parties Alice and Bob 
are able to choose between two observables Ai,A 2 and Bi,B 2 , then one finds that the 
quantity 



/c=0 



{+[P(A 1 = B x + k) + P{B 1 = A 2 + k + l) 
+P{A 2 = B 2 + k) + P{B 2 = A x + k)] 
- [P(A 1 = B x - k - 1) + P(B 1 = A 2 -k) 
+P{A 2 =B 2 -k-l) + P{B 2 =A x -k-l)]} 

with P{A a = (B b + k) mod d) = X]^=o P(A a = (j + k) mod d, B b = j) is bounded by 
the value of 2 for local-realistic theories. In the case of d = 2 this Bell inequality is 
equivalent to the well known CHSH-Bell inequality. 
For the maximally entangled states 



I d-l 

l^o.o) = ^5Zl s >A® \ s )b (H; 
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the parameterization combined with the Nelder-Mead method yields the same values 
up to the order 10~ 5 as found in Ref. [2]. Hence, with increasing dimensionality d, 
the maximum of Id increases from 7 2 = 2.82843 to an upper bound of Id = 2.96981 
for large d. By maximizing the largest eigenvalue of Bj d we obtained that for d > 3 
the inequality is stronger violated by a non-maximally entangled state than by the 
maximally entangled one, which is in compliance with results of Ref. [8]. It should 
be mentioned that a similar result was also found in Ref. [9], where the CHSH-Bell 
inequality was analyzed for neutral entangled K-mesons. These two-state systems 
exhibit a dynamical nonlocality due to their strangeness oscillation and their decay 
property. In particular, it was proven that a non-maximally entangled state violates 
the Bell inequality maximally, whereas the maximally entangled spin singlet state does 
not violate the Bell inequality at all. 

5. The state space geometry of the CGLMP— Bell inequality for bipartite 



In the following we focus on the state space of locally maximally mixed states. As it 
was shown in Refs. [T4"| [15] any locally maximally mixed state can be brought into the 
form 



since any local unitary transformation Ua <8> Ub does not alter the separabil- 
ity/entanglement and (non-)locality properties. Accordingly, these properties depend 
upon only three parameters c,6K which can be regarded as components of a three di- 
mensional vector c. The components of this vector are restricted to lie within a regular 
tetrahedron given by the inequalities 



resulting from the non- negativity condition p > 0. The tetrahedron is spanned by the 
projectors of the four well-known Bell states \^ + ) , |^~) , |$ + ) and |$~). 

For d = 2 it is proven that the PPT (positive under partial transposition) criterion 
is necessary and sufficient, i.e. iff all eigenvalues of a partially transposed density matrix 
are non- negative then the state is separable (for d > 3 it is only a necessary criterion). 
By means of this criterion one finds that all separable states lie within the octahedron 
spanned by the points Ci/ 2 = (0, 0, ±1), C3/4 = (0, ±1,0) and c 5 ~/ 6 = (0,0, ±1) which 
represent equally weighted mixtures of two Bell states. 

Now we want to investigate which locally maximally mixed states violate the 
CHSH-Bell inequality (CGLMP-inequality for d = 2). This inequality is one of the few 



qubits 




(12) 



i=l 



1 + ci + c 2 - c 3 > , 
1 - ci - c 2 - c 3 > , 
1 + ci - c 2 + c 3 > , 
1 - ci + c 2 + c 3 > 
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l*-><*i 



l*+><*+l 





l*-><*-| 



Figure 1. (Color online) Illustration of the state space of all locally maximally mixed 
bipartite qubit states. The connection of all Bell states, located in the corners of the 
cube form a tetrahedron (green). All points inside correspond to states (positivity 
condition). Any equally weighted mixture of two Bell states is separable. Convex 
combinations of these form the octahedron (blue) of separable states. By means of the 
PPT criterion one finds that all states outside this octahedron are entangled. With 
our method we confirmed that states beyond the meshed (red) surfaces violate the 
CHSH-Bell inequality. 



Bell inequalities for which a set of "maximizing" observables can be derived analytically 
[23J . In our case a state expressed by c violates the Bell inequality if and only if at least 
one of the three inequalities 



is satisfied. This means that each of the inequalities defines a cylinder with radius 1 and 
a density matrix that lies outside is nonlocal. The union of the exterior regions of these 
three cylinders belongs to the set of nonlocal states. With the introduced numerical 
method we find up to a high numerical precision the same result. The state space is 
graphically visualized in Fig. [TJ It descriptively demonstrates that not all entangled 
states violate the CHSH inequality. However, it is known that for bipartite qubits all 
entangled states can be distilled and, therefore, contain so-called hidden nonlocality 



c\ + c\ >1 , 
c\ + c\ >1 , 
c\ + c\ >1 



[30j 
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6. The state space geometry for bipartite qudits, the magic simplex W 

For bipartite qudit systems the state space of locally maximally mixed states is more 
complex. This becomes obvious if one tries to find a similar geometric object as the 
octahedron for higher dimensions. Without loss of generality, any locally maximally 
mixed state can be expressed in the form 

d 2 -i \ 

i 3 + ^ m r n ®r m ) (13) 

n,m=l / 

where {Ti, . . . , T (i 2_ 1 } is a set of traceless TrTi = and orthogonal Tr(T\Tj) = kijSij 
operators. As geometric considerations require a real vector space with t nm e R these 
operators also have to be hermitian. The only known candidates for {Ti, . . . , T d 2_ 1 } 
fulfilling all these requirements are the generalized Gell-Mann matrices [31] which 
coincide with the Pauli matrices for d — 2. When they are used to express p one is 
confronted with the problem that in general for d > 2 the expression (fT3|) cannot be 
brought to the form 




d 2 -i 

l> = ^ I 1 3 ® 1 3 + Yl CnTn ® r " I ( 1 1 5 



1 



by means of a certain local unitary transformation Ua®Ub- Furthermore, whereas for 
d = 2 the eigenvectors of p = |(1 <g> 1 + J2n=i c n,o~ n ® o~ n ) are independent of the choice 
of c, for qudits this is not the case. This implies that for d > 2 it is very demanding to 
determine the permitted ranges of the d 2 — 1 quantities c n . 

As one can conclude from this, the geometric structure of the state space of bipartite 
qudits (d > 2) is far more complex than the one of bipartite qubits. Nonetheless, as 
shown in Refs. piJJ [T7J [18] it is possible to construct a simplex spanned by the projectors 
of d 2 orthogonal maximally entangled two-qudit states which can be regarded as a 
generalization of the tetrahedron. This so-called "magic simplex" is given by the set of 
states 

d-l d-1 
k,l=0 k,l=0 

where P^i = \Qk,i) (^k,i\ are the projectors of the d 2 Bell type states produced by 
applying the Weyl operators 

d-1 

w M = Y e¥sfc \ s ) ((* + mod d \ ( 16 ) 

with k, I G {0, d — 1} on the maximally entangled state 

d-1 

\ Q o,o) = ^=Y.\ sA )®\ sB ) 



Vd 

v s=0 
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Figure 2. (Color online) Illustration of the finite discrete classical phase space for 
d = 3 of the locally maximally mixed states of the simplex W. All possible complete 
lines through the point (0, 0) for d = 3 are drawn, representing one class of states 
which have the same geometry concerning separability and (non-)locality. Lines can 
be completed by points with k, I £ {0, .., d — 1} because of the periodicity of the Weyl 
operators implying P k j = Pk+ m -d,i+n-d for all n,m £ Z. 



a 



2 

o 



in the way 

= (W k ,i ® 1) |fi ,o> • (18) 

The simplex is a convex set located in a d 2 — 1 dimensional hyperplane in a d 2 dimensional 
real vector space of hermitian operators spanned by the operators {Pk,i}- It is a subset 
of the set of locally maximally mixed states and the only pure states in W lie at the 
vertices. A detailed discussion on properties of W and on how it is embedded in the 
whole state space can be found in Refs. [161 HZ]- For the subsequent investigations it 
only should be mentioned that local (anti-)unitary transformations Ua®Ub mapping 
W onto itself play an important role for the investigation of separability, entanglement 
and (non-)locality since equivalence classes of states 

[ P ] = y e ny = u A ® u B P u\ ® u j B } 

have the same properties. These equivalences can best be studied when the group 
structure of the Weyl operators and the concept of a finite discrete classical phase space 
are used, see Fig. 2. As shown in Refs. [TBI E2]> states with components {c^} and 
Wk' i'} belong to the same equivalence class if and only if there exists a phase space 



transformation of form 




(19) 



with det = 1 or det = d—1 such that Cki = c k , with m, n, p, q,j,r G Z. 

\p q J \p q J ' 

By means of this concept one can reveal the high symmetry inside W. 
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As there are no operational criteria for bipartite systems with dimension d > 2 that 
can solve the separability problem completely, one uses different methods such as the 
PPT criterion, constructing optimal entanglement witnesses or applying entanglement 
measures. In precedent publications on the geometry of separability of W [HI [T71 [T51 S3] 
these methods were combined to determine regions of separable and entangled states. 
Some known results for d = 3 will be used in the subsequent discussion. 

In this present work we want to determine the set of states p of the magic simplex 
W that violate the CGLMP-Bell inequality by using the previously introduced method. 
In particular we are interested in the set of states which lie on the boundary, i.e. for 
which 

max Tr(B Id p) = 2 

holds with high numerical precision. Determining states with this property is in gen- 
eral computationally intensive. In Ref. [32] the following theorem was stated and proven: 



Theorem 1: If for a given state r the maximal value of maxi^r) = max^ Tr(Bi d r) 
is known then the state p = ^jr-ld 2 + v T with added uncolored noise ^r\& lies on the 
boundary of CGLMP-Bell violation for v = 2/ max/j(r). Proof: The proof is based 
on the fact that the trace of the CGLMP-Bell operator vanishes, i.e. Tr(Bj d ) = which 
leads to 



max Tr(Bi d p) = max Tr{B 



1 - V 



) 



max Tr{Bi d v r) 



2 <^> v = 2/max/ rf (r) . 



7. The state space geometry of the CGLMP-Bell inequality for bipartite 
qutrits 

In the following we investigate and illustrate low- dimensional sections of the magic 
simplex for bipartite qutrits (d = 3). For the numerical optimization of the Bell operator 
the Nelder-Mead method [2S] was used, which performed relatively fast and reliable. As 
a first example consider the so-called isotropic states 

p=^^l 9 + aP k , l . (20) 

The phase space transformation rules (TTOT) imply that all such one parameter states with 
arbitrary k, I G {0, ..,2} but same a have the same properties in terms of separability, 
entanglement and (non-) locality. For this set the optimal entanglement witness yields 
that p is separable for | > a > and entangled for a > |. Numerical investigations are 
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in compliance with the conjecture made in Ref. [2] that the CGLMP-Bell inequality is 
violated for a > |(6v^3 — 9). 

Next, consider the two-parameter families of states of the form 

p= 1 ~ (X ~^ tg + aP k>l + f3P m , n (21) 
y 



with arbitrary k,l,m,n G {0, ..,2}. Any such state is local unitary equivalent to the 

-Q- 

"IT 

8a 2 + 8(3 2 - llfia + 2a + 2(3 - 1 = . 



state p = 1 " - lg + aPofl + /3Po,i- For this set of states the PPT boundary derives to 



There are also bound entangled states in this set. They can be found through optimal 
entanglement witnesses yielding the boundaries [33] 

4a 2 - 5a + 40/3 2 + (17a - U)(3 + 1=0 and 
A(3 2 -5(3 + 40a 2 + (17(3 - 14)a + 1 = 0. 

The boundary of the CGLMP-Bell violation was obtained by computing max Iz(p) for 60 
equally spaced points on the boundaries of positivity (a = ^p, (3 = 9 -^- and (3 = 1 — a). 
Theorem 1 was then exploited to determine the values of a and [3 corresponding to states 
on the boundary, i.e. max^ Tr {Bi d p) = 2. The state space is graphically illustrated 
in Fig. [3J 

The illustration suggests that the boundary of CGLMP-Bell inequality violation 
describes a circle for a, (3 > and a line if one of the parameters is negative, i.e. a < 
or (3 < 0. Note that this abrupt change in the shape only appears for the boundary of 
CGLMP-Bell inequality violation but not for the boundary of separability. 

Conjectures on the exact specifications of these boundaries are contained in the 
considerations of the following three-parameter families of states 

1 — a — 6 — t 

P = ^ 1 1 + aP Kl + /3P m , n + 7 P M (22) 

y 

as the special case 7 = 0. Depending on whether the index pairs {(k,l), (m,n), (p,q)} 
form a line or not this three-parameter family of states is either local-unitarily equivalent 
to the state 

1 — a — (3 — 7 

p line (a, (3, 7) = ^— '-1 9 + aP ,o + 0P o ,i + 7^o,2 (23) 

or 

Poffiinc(a,/3,7) = 1 ~ a ~ l 9 + aPo,o + PP ,i + 7^i,o • (24) 

y 

The class of line states are graphically illustrated in Fig. H] and a particular slice of the 
off line states is visualized in Fig. [5] (b). 
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bound entanglement 



Figure 3. (Color online) Illustration of the state p = 1 9 + aP 0t0 + PPo,i- All 

states lie within the (green) triangle which corresponds to the border of positivity. 
The (blue) ellipse corresponds to the PPT border, i.e. all PPT states lie inside. As 
shown in Ref. [IT if one parameter is negative, there is also a small (grey) filled region 
of bound entanglement. States p in the (red) filled area violate the CGLMP-Bell 
inequality (max Is(p) > 2). Interestingly, the geometry given by the CGLMP-Bell 
operator changes its shape (from a circle to a line) at the transition from positive to 
negative parameters. 

For the states pn ne ((x, ft, j) the PPT boundary is 

8a 2 + 8/3 2 + 87 2 + 2a + 2(3 + 2 7 - lla/3 - lltry - 11/37 -1=0. 

The separablity boundary obtained by optimal entanglement witnesses revealing further 
(bound) entangled states is 

40a 2 + a(17/3 + 177 - 14) + 4/3 2 + 7 (4 7 - 5) - 0(197 + 5) + 1 = 

as well as all possible permutations with respect to a, /?, 7. The boundary of CGLMP- 
Bell violation is deduced from the values max/ 3 (p) of 920 equally spaced states on the 
boundaries of positivity (a = ^"g" 1 (and all parameter permutations of this term) and 
7 = 1 — a — /3). The resulting points apparently describe a spherical surface which meets 
the boundaries of positivity and PPT at the point a = /3 = 7 = ~ and is intersected by 
planes in the vicinity of the isotropic states. 

When taking into account symmetries and the fact that the suggested analytical 
ideal measurements in Ref. [2] yield the boundary parameter |(6a/3 — 9) for the isotropic 
states one can suppose a sphere with radius r = yfg(413"\/3 — 558) and central point 
having the coordinates a = (3 = 7 = ^(—361 + 186\/3). These specifications coincide 
with the numerical data up to the order 10 -5 . Presumably, these discrepancies decrease 
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Figure 4. (Color online) Illustration of a three dimensional state subspace given 
by bipartite qutrit line states pu ne (a, (3, 7), Eq. (f23|) . As for bipartite qubit states, 
Fig. [U all states have to be within a (green) tetrahedron (positivity condition). The 
boundary of PPT states forms a (blue) cone and thus all states beyond this surface 
are entangled. The tip of the cone touches the surface of positivity ata = /3 = 7= |. 
States p beyond the shaded (light red/blue) surface violate the CGLMP-Bell inequality 
(max 23 (p) > 2). Again, the geometric form changes from a sphere to a plane when 
one or more parameters become negative. 



when the numerical optimization method is set to a higher working accuracy. The planes 
are given by the function 

7= ^(a + (3 + 6V3 - 9) (25) 

and the two additional functions given by the permutations with respect to a, (3 and 7. 
In contrast to the specifications of the sphere, these functions are not conjectured but 
result from the measurement settings given in Ref. [2] and local unitary transformations 
constituted by the Weyl operators. 

For the remaining family of the off line states, Eq. f[2~4"|) . the geometric form of the 
boundary of CGLMP-Bell violation has a more complex shape. An uncovered suggestion 
on the exact form has therefore not been made. Fig. [5] (b) gives an impression of the 
fact that the boundaries of PPT, as well as CGLMP-Bell violation have a different 
shape for the different types of states. For instance, in contrast to the mixtures of states 
on a line where there is only one state a = /3 = 7 = |on the boundary of positivity 
1 — a — (3 — 7 = that does not violate the CGLMP-Bell inequality, for the second 
type there is a whole region of states for 1 — a — (3 — 7 = that obeys the CGLMP-Bell 
inequality. Moreover, in comparison the entire region of nonlocal states in Fig. |5] (b) is 




Figure 5. (Color online) Comparison of the class of line states pn ne (a, lr, f )' Eq. (|23|) . 
with the class of off line states /0 o ff Unc(a, §, §), Eq. (l2~31) . In the first case, one finds 
that there is no bound entanglement and that the separable states form a polygon. In 
the second case, the boundary given by the PPT criterion is curved (blue) and bound 
entanglement is found for (3 < (grey filled region). The region of nonlocal states 
(filled red region) is smaller in case the third Bell state does not belong to the line 
formed by the other two Bell states of the mixture. In Sec. [8] an entanglement measure 
for this subset of line states is derived, for which it turns out that for p\i nc (a, §, §), 
all states on lines parallel to the border of separability contain the same amount of 
entanglement. This proves that the amount of entanglement is not related to the 
violation of the CGLMP-Bell inequality. 



smaller than in Fig. (a). For entanglement, one finds the contrary, which is a rather 
counterintuitive result. This is due to the fact that the region of separable states of the 
subset of line states 

PlineK ^ , ^ ) ( 26 ) 

already is larger than the PPT region of the same subset of off line states 

Poffim e (a,-,-) (27) 
which in addition also contains bound entangled states for (3 < 0. 



8. The amount of entanglement compared to Bell— inequality violation 

In Ref. [19J a general entanglement measure for multipartite systems of any dimension 
was introduced as well as a method how bounds can be derived for it. Here we apply 
this measure to bipartite qutrits. 



On the state space geometry of the CGLMP-Bell inequality 



16 



For any bipartite qutrit state p = "^iP^ipi) the entanglement measure is given 

by 

£{p):= inf Y,Pi {S(Tt a \^)(^\) + S(Tr B \^)m} 

I 

= 2 inf V^5(Tr A |^)(^|) , 

I 

where S is any Renyi entropy and in this case we choose S = — log 2 (l — 8/9 Si), where 
Sl is the normalized linear entropy, i.e. Sl(p) = f (1 — Trp 2 ). In general the infimum, 
i.e. the optimal decomposition is not known. The derivation of lower bounds for the 
measure are based on the observation that the entropy can be re-written by an operator 
sum, where the single terms are called general concurrences or m-concurrences C as 
they are measures in their own right. For bipartite qutrits we obtain 

S{p A ) = S(p B ) = -21og 2 (l - \{C\p A ))) , (28) 

where p A and pB are the reduced density matrices. The concurrences C are for any pure 
state ipi given by 



1/2 

(29) 



For any mixed state p the following bound can be obtained 



C 2 (p) > ~5>ax jo, (2max{A*} - £ A*) J 



2 



(30) 



where A^ are the squared root eigenvalues of 

pO k p*(O k y. (31) 



In the case of bipartite qutrits via Eq. fl28|) one obtains nine operators O k , which can 
be chosen to be the tensor products of all antisymmetric Gell-Mann matrices. Clearly, 
entanglement does not change under local unitary operations U a ®UbpU\® U^, but 
this does in general not hold for the bounds. In some cases the bounds can be further 
optimized by two elements U A and Ug of the SU(d) group (d = 3 in this case). Here, 
too, one can use the generalized Euler angle parameterization. 

Consider the class of line states where two Bell states are equally weighted, i.e. 
Eq. ( 126|) . For these states bounds on the entanglement measure were derived in Ref. [19J. 
It was found that for a > \ + | the bounds agree with the PPT criterion and as there is 
no bound entanglement in this class of states Ref. [T7] this means the measure is exact. 
However, for a < \ + f the bounds were found to be worse than the border given by 
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PPT. By optimizing Ua and Ub with the proposed method we obtain explicitly: 



1 P 
For a > - + — 

~ 4 8 



S(pi inc (a 




' 2' 2 



)) =-21og 2 |l- I ^max{0,-(2-8a + /3)} 2 | 



1 P 
For «<- + —: 

4 8 



(32) 



(33) 



£ OineK f , §)) = -2 log 2 |l - ^ max{0, -(2 + 4a - 5/3)} 2 | 

which exactly coincides with separability. In this sense, the entanglement measure is 
exact for the whole class. 

It implies that all states lying on a line parallel to the separability boundary have 
the same amount of entanglement. However, although the amount of entanglement is 
the same for these states, the maximal achievable value of I3 varies significantly. As 
an example consider the family of states pn ne (a, 8a ~ 5 , 8o T 5 ) for which, according to 
Eq. (1321) . the entanglement measure yields the same value. Depending on a 6 til' 1] 
one finds the following values for maxl^: 



a 


0.611 


0.622 


0.633 


0.644 


0.656 


0.667 


max/3 
viol. [%] 


1.915 



1.819 



1.875 



1.993 



2.123 
6.170 


2.265 
13.226 



This reveals that there exists no simple connection between the amount entanglement 
and the strength of CGLMP-Bell inequality violation. Another example is illustrated 
in Fig. where the amount of entanglement of states on the border of CGLMP- 
Bell inequality violation (max Is(p) = 2) and of positivity are plotted. Both results 
demonstrate that there exists no monotone function relating the Bell inequality violation 
to the amount of entanglement. 



9. Summary 

We presented a detailed analysis of the CGLMP-Bell inequality in terms of state 
space geometry for qubits and qutrits. It was necessary to ascertain whether this Bell 
inequality is violated or obeyed for a given quantum state p. In order to solve this, 
we used a reduced form of the Euler angle parameterization of the SU(d) group and 
numerical methods. By means of this parameterization introduced in Sec. [3] we were 
able to avoid problems with constraints and to reduce the numerical search to d 2 — d 
real parameters for each observable used in the Bell inequality. In general, this method 
can be adapted to all kinds of qudit systems and, moreover, has the advantage of not 
being restricted to a certain Bell inequality. 

In Sec. H]we tested the proposed method for the CGLMP-inequality and compared 
it with some known results. We found that the method works within good numerical 
precision. 



On the state space geometry of the CGLMP-Bell inequality 18 




0.2 0.4 0.6 0.8 1.0 



Figure 6. (Color online) These curves show the amount of entanglement for the line 
state p(a, §, §), Eq. in dependence of j3 (compare with Fig. [5] (a)). The solid 

thin (black) curve is the amount of entanglement for all states lying on the boundary 
of positivity. The other graphs show the amount of entanglement for the boundary 
of CGLMP-Bell inequality violation. In detail, the dotted (red) graph illustrates the 
amount of entanglement for /3 < 0. The solid thick (green) graph illustrates the amount 
of entanglement in the range j3 = to the separable state p(h, |, |) along the curved 
CGLMP boundary, while the dashed (blue) graph illustrates the region /3 > |. 

In Sec. [5] we began our geometrical state space considerations concerning the 
CGLMP-inequality. Here we investigated the tetrahedron spanned by the four well- 
known Bell states which represents the state space of locally maximally mixed states for 
bipartite qubits. We determined the sets of entangled and CGLMP-inequality violating 
states with analytic methods and the numerical method again yielded the same result. 
We visualized the results in Fig. [TJ The illustration descriptively shows that not all 
entangled states violate the CGLMP-inquality. 

In Sec. [6]we looked at the state space geometry for the more general case of bipartite 
qudits. We focused on states which are elements of the so-called "magic" simplex W - 
a generalization of the tetrahedron. The simplex W consists of all convex combinations 
of d 2 maximally entangled bipartite qudit states which can be obtained by choosing one 
maximally entangled state and then applying the Weyl operators in one subsystem. All 
states of the simplex have the property of being locally maximally mixed, i.e. any partial 
trace results in a maximally mixed state. By group theoretical methods certain classes 
of states within W can be found which have the same geometry concerning separability, 
entanglement and (non-)locality, see Fig. [TJ 

In Sec. [7] we focused on revealing the nonlocal states of the "magic" simplex W for 
the special case of bipartite qutrits (d = 3). We exploited a particular property of states 
containing uncolored noise (1^) to further reduce the complexity of the optimization 
(Theorem 1). We found a large region of entangled states which does not violate 
the inequality. In particular, the boundaries of CGLMP-Bell inequality violation and 
separability are visualized for the class of the so-called "line" states and for the class 
of the so-called "off line" states in Fig. [3] to Fig. [5j As could be seen in all cases, both 
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boundaries are not simply related by a constant shift or scaling factor but rather are 
fundamentally different. 

In Sec. M we considered the relation between a recently introduced (multipartite) 
entanglement measure p35] and Bell inequality violation. For this measure, we derived 
lower bounds for the class of line states with two equally weighted parameters. Via 
the proposed optimization method we were able to optimize these bounds. Comparing 
the strength of CGLMP-Bell inequality violation and the amount of entanglement, 
we explicitly demonstrated that in general there is no monotone relation between the 
amount of entanglement and a Bell inequality violation. This extends the result found 
in Ref. [8] from pure to mixed states and strengthens the conjecture that nonlocality 
and entanglement are different features. 
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